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We show how recent state-reconstruction techniques can be 
used to determine the Hamiltonian of an optical device that 
evolves the quantum state of radiation. A simple experimen- 
tal setup is proposed for measuring the Liouvillian of phase- 
insensitive devices. The feasibility of the method with current 
technology is demonstrated on the basis of Monte Carlo sim- 
ulated experiments. 

PACS: 03.65.Bz, 43.55.-i 



In recent years, the possibility of "measuring" the 
quantum state, after remaining for a longtime a mere 
theoretical speculation eventually entered the realm 
of true experiments. From the first experimental demon- 
stration Q, the so called "homodyne tomography" tech- 
nique advanced to the level of a quantitative state- 
reconstruction technique BEj , achieving a high degree of 
reliability in experiments g|. This state- reconstruction 
method is now ready to be used for concrete applications. 

What is the practical use of measuring a quantum 
state? Apart from the availability of a kind of "univer- 
sal detector" |j) that provides informations on all ob- 
servables at a time, measuring a quantum state is the 
only way to check a state preparation within a (gener- 
ally nonorthogonal) set. In turn, the use of homodyne 
tomography to test state-preparation becomes a way to 
check the operation of a quantum device that prepares 
a chosen state from a given one. It is now natural to 
ask if eventually it would be possible to recover a com- 
plete information on the quantum device itself, namely 
to reconstruct the detailed form of its Hamiltonian — or, 
more generally, of its Liouvillian, as in reality the device 
is always an open quantum system. Previous theoretical 
proposals to give a complete characterization of quan- 
tum processes have been made in Refs. j^). There, the 
methods are restricted to systems with finite dimensional 
Hilbert space, and the method does not lead to an ex- 
plicit reconstruction of the Liouvillian. In this letter we 
show how this goal can be achieved in practice, presenting 
a simple experimental setup for measuring the Liouvil- 
lian of a phase-insensitive optical device, using currently 
available technology. 

The main idea for reconstructing the Liouvillian of a 
quantum device is sketched in Fig. |l|. One should im- 
pinge the device with a known input state pi n from a 
(over) complete set, then determine the state p ou t at the 
output, and finally compare pi„ to p ut- For an opti- 
cal device the determination of the output state is made 
possible by the homodyne-tomography technique. Re- 
garding the generation of the set of input states {p, n }, 



an experimental method is suggested later in this letter. 
The evolution of the state from pi„ to p ou t is governed 
by the Green (super)operator Q 



Pout — Qpi 



(1) 



where Q has actually a four-index matrix representa- 
tion, and on the Fock basis one has (n\p ou t\m) = 
£~ fc=0 G*t{h\ Pin \k). 
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FIG. 1. Sketch of the method for measuring the Liouvillian 
of an optical device. A known input state p;„ is impinged into 
the device, and a quantum tomography of the output state 
is performed using a homodyne detector. By scanning an 
(over) complete set of states pi„ at the input and comparing 
them with their respective output states, it is possible to re- 
construct the Liouvillian of the device. The histograms of 
homodyne data, here given for the sake of illustration, corre- 
spond to a device that consists of an empty cavity, and with 
the input states as number states pi„ = \n){n\. 

For a device that is homogeneous along the direction of 
light propagation the Green superoperator can be written 
as the exponential of a constant Liouville superoperator 
C as follows 



Q = cxp(£r) 



(2) 



where r is the propagation time (i.e. the device length). 
The Liouvillian C gives the evolution of the state through 
an infinitesimal slab of the device media according to 
the master equation p = Hp. In this letter we restrict 
our attention to the case of a perfectly phase-insensitive 
device: as it will be clear from the following, the case of 
phase-sensitive device is much more complicate, and will 
be analyzed elsewhere [||. A phase-insensitive device, is 
a device that leaves dephased states as dephased, as in 
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the case of a traveling wave laser amplifier. A dephased 
state is diagonal in the photon-number representation, 



with density matrix of the form p = J2r 



\n)(n\, 



where {\n}} denotes the complete set of eigenvectors of 
the photon- number operator a^a of the field mode with 
annihilation operator a. For the evolution of dephased 
states it is sufficient to determine the sector of the Green 
superoperator that evolves dephased states, i.e. the two- 
index Fock matrix G nm = (n\ Q[\m)(m\] \n). 

The experimental reconstruction of L nm could be per- 
formed by impinging a number state p in — \n)(n\ on the 
device, and then making the homodyne tomography of 
the output state. In this fashion, the number probabil- 
ity distribution rk(n) of the output coincides with the 
n-th row of the Green matrix G n m, and by varying n 
one would reconstruct the whole matrix. Since produc- 
ing number states is experimentally difficult, one would 
try to use coherent states instead. In this way matrix el- 
ements of the form Q[ \a)(a\ ] \ip') would be obtained, 
with \a) denoting the scanning coherent input state, 
and \ip) and being a couple of vectors of the tomo- 
graphically reconstructed matrix representation. Unfor- 
tunately, the relation between G nm and (n| Q[ \ct)(a\ ] \n) 
is highly singular, involving the P-function of |m)(m|, 
and hence the matrix G nm cannot be obtained in this 
way starting from experimental data. On the other hand, 
the Fock representation has a privileged role, because 
here the Liouvillian matrix has a transparent meaning 
in terms of creation and annihilation operators. How 
to overcome the problem of generating input number- 
states? Actually, for our purpose, it is sufficient to gen- 
erate number states with random n as far as n is known. 
This leads us to devise the setup depicted in Fig. 0. 
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FIG. 2. Experimental setup, including the apparatus used 
to generate the input number states p in = |n)(n needed for 
the tomographical reconstruction of the Liouvillian of an op- 
tical device. A random-n Fock state \n) for the input beam 
is achieved by performing photodetection at D on the other 
twin beam, n being the measured number of photons. A non 
degenerate optical parametric amplifier (NOPA) with vacuum 
input is used to produce the twin beams. 

A non degenerate optical parametric amplifier (NOPA) 
with a strong classical pump down-converts the vacuum 
into a pair of twin beams. The twin beams are used 
as a random-n Fock state generator, by measuring the 



number of photons on one beam (detector D in Fig. 
||) while impinging the other beam on the optical de- 
vice. For quantum efficiency T)d = 1 at detector D, 
the photodetection would reduce the twin-beam state 
|TB) cx X)n=o n ) m to a random-n state \n) at the 
input of the optical device, with thermal probability dis- 
tribution w n — \k\ 2ti (l — |k| 2 ), where n is the measure- 
ment outcome at D. The tomographically reconstructed 
number probability (k\ Q[ \n)(n\ ] \k) of the output state 
already would provide the nth row Gkn of the Green ma- 
trix. On the other hand, for rjp < 1, a mixed state p n 
will actually enter the device instead of \n)(n\, as a re- 
sult of state reduction at D. The outcome n probability 
distribution then becomes 



Pn = (1 - M 2 ) 



2\n 



[( Vd -1)\k\ 2 + !]»+! 



(3) 



One can easily show that the tomographically recon- 
structed output number probability rfe(n) = {k\Q[p n \ \k) 
is related to the Green matrix through the identity 



r k (n) = [( VD - 1)\k\ 2 + l] n+1 
x { m + n 



m=0 



[\K\ 2 (l-r, D )} m G k , m+n 



(4) 



The relation (W) can be inverted as follows 



G k 



ki 



[(^-i)M 2 + i] 
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n + I 
I 



(rj D - 1)\k\ 2 + 1 



r k (n + l). (5) 



Eq. (||) is our algorithm to reconstruct the Green matrix 
Gki from the collection of all tomographically measured 
number probabilities rfc(n) for different outcomes n (in 
practice the sum in Eq. (|5|) is truncated at some max- 
imum n). Notice the interplay of the gain k and the 
quantum efficiency j]d in determining the probability p n 
on one hand, and in producing the statistical errors in 
the reconstructed Green matrix on the other hand. For 
decreasing quantum efficiency rjr> — ► 0, larger values of n 
can be made more probable by increasing the gain of the 
NOPA as k — > I~. However, at the same time, conver- 
gence of the sum in Eq. (H) becomes slower, and statis- 
tical errors of matrix elements Gki increase as result of 
tomographic errors on r k {n). Hence, the effect of quan- 
tum efficiency rju, which reduces the size of the viewable 
matrix Gki, can be partially compensated by increasing 
the gain of the NOPA, however at expense of statistical 
errors for Gki- For the tomographic measurement, by in- 
creasing the number of experimental data and using T]h- 
dependent pattern functions, the method of Ref. ^ can 
compensate the effect of low quantum efficiency r\u < 1 , 
that, anyhow, must be above the threshold i]h = 1/2. 
On the other hand, for quantum efficiency r\B at detec- 
tor D there is not such a threshold, as one can see from 
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convergence and and error-propagation analysis of Eq. 

The proposed state-reduction scheme — based on twin- 
beams from a NOPA — is not a new one, and, for example, 
a similar setup has been proposed in Ref. (To) to gener- 
ate Schrodinger-cat states. As such state-reduction is the 
core of our measurement method, we want to examine it 
at work in a realistic situation. Typically the NOPA 
can be pumped by the second harmonic of a Q-switched 
mode-locked Nd:YAG laser, with the output twin-beams 
pulsed at a repetition rate of 80 MHz, and with a 7ps 
pulse duration. Thus, the twin-beam mode with anni- 
hilator a at the input of the optical device is actually a 
wideband mode, with frequency centered around 532nm, 
and width of 140 GHz (the inverse of the pulse time- 
length). The same Nd:YAG laser beam is used for the 
local oscillator (LO) of the homodyne detector H. In 
this way the LO has the same central frequency and the 
same time-envelope of the beam at the input of the op- 
tical device. The integration time at photodetector D 
can be set to 1 ns, which is greater than the pulse width 
and shorter than the distance between pulses. In this 
way each pulse is completely annihilated by the detector 
D during the integration time, and, correspondingly, the 
homodyne measurement is made with the LO matched 
on the same pulse shape of the signal twin-beam, which 
means that the measurement is performed on the right 
wideband mode after reduction by detector D. Moreover, 
the detector D and the optical device can be geometri- 
cally placed in such a way that, within a narrow solid an- 
gle, the direction of the respective input fc-vectors are the 
same relative to the k- vector of the NOPA pump, so that 
state reduction at D affects only radiation at the twin 
k at the input of the device, so that the state-reduced 
modes at D and at the device are perfectly matched. 

From the above scenario it follows that the mode with 
annihilator a at the input of the optical device is actually 
a wideband mode, and hence we measure the effective 
Liouvillian over a 140 GHz bandwidth centered around 
532nm. Then, it is clear that all measurements for differ- 
ent random inputs can be considered as independent only 
if the (atomic) relaxation times in the optical devices are 
shorter than the pulse-repetition period. 

Now we show results from some Monte-Carlo simulated 
experiments to see our method at work, and estimate the 
number of measurements needed for the reconstruction of 
Gki- 

The simplest phase insensitive device is the phase- 
insensitive linear amplifier (PIA), with Liouvillian C 



C = 2{AV[a)] +BV[a}} 



(6) 



where V[6]p = 6p0^ - l/2(0t0p + pO^O) for any complex 
operator 9. The Liouvillian matrix has the form 



(n\C 



\m){m\ \n) 



(7) 



where So- denotes the Kroneker delta. Notice that L nm is 
tridiagonal, the upper diagonal corresponding to the one- 
photon absorption apa^ of the loss term V[a], the lower 
diagonal corresponding to the one-photon emission pa 
of the gain term 2? [at], and the main diagonal containing 
the anticommutators coming from both terms. 
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FIG. 3. (a) Linear phase insensitive laser amplifier with 
A = .1940 and B = .00945. Monte Carlo simulation of the 
reconstruction of the Liouvillian L nm in Eq. |^ by means of 
the proposed experimental setup. The reconstruction is per- 
formed by using 2 statistical blocks of 10 6 homodyne data for 
each of the output states (from a total of 2.7 * 10 11 of data). 
Quantum efficiencies are 770 = .8 and r\u = .85. The gain of 
the NOPA is K = .6. (b) the same reconstruction for lower 
quantum efficiency 770 = .3. Here the gain is set to n — .4, 
and the same number of data is used out of a total of 9.9 + 10 11 . 




2 jj4(m + l)[6 nm+ i - S nm ] + Bm[5 nm -i - £„,„]} 



FIG. 4. The three main diagonals of the Liouvillian L„ m 
of Fig. ^|a are given with their statistical error bars. The full 
line is the theoretical value from Eq. (Q). 

In Fig |]a we show a typical result of a Monte Carlo 
experiment of the tomographic reconstruction of the Li- 
ouvillian (Q). We used quantum efficiency tjd — .8 and 
rjD = -3 at detector D and r\n — .85 at the homodyne de- 
tector H. One can see that the details of the matrix are 
well recovered, and the truncation of the Hilbert space di- 
mension does not affect the reconstruction. In Fig. |^ the 
three main diagonals of the matrix are plotted with their 
statistical errors against the theoretical value, showing a 
very good agreement. The statistical errors are of the 
same size of those of the tomographically reconstructed 
output probabilities. Notice that the number of data 
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10 11 ~r~ 10 12 needed for this experiment could be collected 
in a few minutes at the repetition rate of 80 MHz. In Fig. 
||b we present a simulated experiment using the very re- 
alistic value rjD = -3 of quantum efficiency, however for 
the reconstruction of a smaller matrix 5x5. Notice that 
only the photodetector D is required to be linear single- 
photon resolving, whereas the homodyne detector takes 
advantage of amplification from the LO, and hence can 
use high efficiency detectors (the value r\u — .85 here 
used has been widely surpassed in the real tomographic 
experiments, as in Ref. ||). 

As another example, we simulated the experimental 
tomography of the effective Liouvillian of a one-atom 
traveling- wave laser amplifier. 
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FIG. 5. Theoretical Liouvillian for a one-atom laser, ob- 
tained by solving numerically the master equation (^). The 
parameters for this laser are C = = 12; n a = -^p- = 7; 
a = 1; / = = 1; 7 = 1; U = .0115. 



where g is the electrical-dipole coupling, 7|| and 7j_ are 
the decay rates of population inversion and atomic po- 
larization respectively, 7 is the cavity decay rate, Co is 
the unsaturated inversion (—1 < cr® < 1), a± z are the 
Pauli matrices (with 0, ±1 entries), and p now denotes 
the joint atom-radiation density matrix. In the quan- 
tum jump simulation, the atom is traced out at a time 
t* 7|j~j_- hi Fig- H a Monte Carlo simulated tomo- 
graphic experiment is shown for reconstructing the Liou- 
villian in Fig. [| (the output homodyne probabilities are 
simulated starting from the quantum jump Green ma- 
trix). One can see how the method allows a detailed re- 
construction of L nm , including not only one-photon pro- 
cesses on the three main diagonals, but also multiphoton- 
absorptions on the upper triangular part. 

In conclusion, we have seen that it is possible to ex- 
perimentally reconstruct the Liouvillian of a quantum 
optical phase-insensitive device, using homodyne tomog- 
raphy in a scheme based on parametric down conversion 
from a NOPA. We have shown the feasibility of the recon- 
struction with an experimental setup that uses standard 
technology devices. The problem of low efficiency at the 
single-photon resolving detector D — the major obstacle 
for the experiment — has been solved by implementing a 
compensation algorithm that makes the reconstruction of 
5x5 Liouville matrix possible even for r\D = .3, and with 
a number of data that can be collected in a few minutes 
of experimental run. 




FIG. 6. Monte Carlo simulated experiment for the recon- 
struction of the laser theoretical Liouvillian in Fig. ^. Here 
r]D ~ .8, r\H = .85, and k = .65. A set of 8 * 10 7 homodyne 
data have been used of a total of 3.7* 10 11 measurements with 
random photon number at detector D. 

In Fig. |5| the theoretical Liouvillian matrix is plotted, 
as obtained from a long-run quantum jump simulation of 
the one-atom-laser master equation fll| : 

P = + <ro)T>[a+] + f (1 - ir )V[v-] + 

1 (7-L - 1) »W + lT>[a]\p + g[a+a - a_a), p) (8) 
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